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Asymptotic Satellites near the Straight-Line Equilibrium 
Points in the Problem of Three Bodies. 

By Daniel Buchanan. 



§ 1. Oscillating Satellites. 

It was shown by Lagrange in a prize memoir * in 1772 that if two finite 
spheres revolve in circles about their common centre of mass, then there are 
three points on the line joining their centres at which an infinitesimal body 
would remain if it were given an initial projection so as to be instantaneously 
fixed with respect to the moving bodies. These points are called the straight- 
line equilibrium points of the problem of three bodies. Starting from minus 
infinity the order of the equilibrium points and the finite bodies is an 
equilibrium point, a finite body, a second equilibrium point, the other finite 
body, the third equilibrium point. 

If the infinitesimal body is given an initial displacement from a point of 
equilibrium and initial conditions are so chosen that it moves in a closed orbit 
relatively to the moving system, it is then called an oscillating satellite. 

The problem of the oscillating satellite has been discussed extensively. In 
the papers cited below,t the differential equations are limited to their linear 
terms and the orbits restricted to the plane of motion of the finite bodies. 

A rigorous demonstration for the existence of periodic orbits for the 
oscillating satellite and a practical method for constructing them are given by 
Moulton in Chapter V of his "Periodic Orbits." J In this memoir the differ- 
ential equations are unrestricted in the number of terms which may be taken, 

♦Lagrange, "Collected Works," Vol. VI, pp. 229-324. 

fThe following references to the literature of the oscillating satellite are taken from Chapter V 
of Moulton's "Periodic Orbits": 

Poincare, Les Methodes Nouvelles de la Mdcanique Celeste, Vol. I (1892), p. 159. 

Burrau, Astronomisohe Nachrichten, Nos. 3230, 3251 (1894) . 

Perchot and Mascart, Bulletin Astronomique, Vol. XII (1895), p. 329. Moulton adds, "apparently 
their work is vitiated by an error in establishing the existence of the solutions, and their construction 
fails where they stopped." 

Sir George H. Darwin, Acta Mathematica, Vol. XXI (1897), p. 99. 

Plummer, Monthly Notices, Royal Astronomical Society, Vol. LXIII (1903), p. 436, and Vol. LXIV 
(1903), p. 98. 

$This memoir will be cited as the "Oscillating Satellite." Another method is given in Chapter 
VI of the " Periodic Orbits," but we are concerned with the former method only. 
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and the orbits are not limited to two dimensions as in the previous literature. 
Three classes of orbits are shown to exist and they are designated according to 
their periods as Class A, Class B, and Class C. The orbits of Class A and 
Class C are of three dimensions, while those of Class B are of two dimensions. 
The orbits of Class C are shown to exist under special conditions which might 
never be realized in the problem. Their period is a multiple of the periods of 
Class A and Class B when these latter periods are commensurable. Practical 
constructions are made for the orbits of Class A and of Class B, but, owing 
to the complexity of the problem, no attempt has been made to determine 
whether orbits of Class C exist which are distinct from those of Class 
A and Class B. 

§ 2. Asymptotic Satellites. 

The object of this paper is to obtain solutions of the differential equations 
of motion of the infinitesimal body which will approach the periodic solutions 
of Class A and Class B as the time approaches plus infinity or minus 
infinity. When the infinitesimal body moves in an orbit defined by such 
solutions, it will be called an asymptotic satellite. 

The question of the existence of solutions which are asymptotic to the 
orbits of Class C is not considered in this paper, owing to the fact that it has 
not been determined in the " Oscillating Satellite" that orbits of Class C exist 
which are distinct from those of Class A and Class B. 

The orbits which are asymptotic to the equilibrium points themselves have 
been determined by Warren.* These orbits are of two dimensions and lie in 
the plane of motion of the finite bodies. 

The form of the asymptotic solutions is that adopted by Poincare, f each 
term being of the type e Xt P(t), where ^ is a constant having its real part 
different from zero, and P(t) is a constant or periodic function of t. It has 
been shown by Poincare that solutions of this type will converge for all values 
of t, provided that certain divisors which appear in the construction of such 
solutions are different from zero.$ If these divisors vanish, terms of the form 
te Xt P(t) will arise. If, therefore, the construction can be made so that no 
terms occur in t explicitly, the divisors previously mentioned are different from 
zero and the solutions will converge for all values of t. Hence it is sufficient 

* Warren, " A Class of Asymptotic Orbits in the Problem of Three Bodies," American Joubnai, of 
Mathematics, Vol. XXXVIII, No. 3 (1916), pp. 221-248. 
t Poincare, "M<5canique Celeste," Vol. I, p. 340. 
$ Poincare^ loc. oit. p. 341. 
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to consider only the formal construction of the asymptotic solutions, and if 
solutions can be constructed so as to contain no terms in t explicitly, their 
convergence is assured by Poincare's theorem. 

§3. The Diffebential Equations of Motion. 

Let the motion of the infinitesimal body be referred to a set of rotating 
rectangular coordinates £>?£ of which the origin is at the centre of mass of the 
finite bodies, the £-axis is the line joining the finite bodies, and the ^-plane is 
the plane of their motion. The £- and >7-axes rotate about the £-axis in the 
direction of the motion of the finite bodies and with the same angular velocity. 
The units of length, mass, and time will be taken so that the distance between 
the finite bodies, the sum of their masses, and the Gaussian constant respectively 
shall each be unity. With the units thus chosen, the mean angular motion of 
the system is likewise unity. Let the masses of the finite bodies be denoted by 
1 — (i and (i, 0<^<|. On denoting the coordinates of the infinitesimal body 
by £, v, £> an d differentiation with respect to t by primes, the differential 
equations of motion are * 

*. o*'- dU y,"+2P- dU Z"- dU 



2Z7=f+) 7 2 + 2(1 ~' f) + — 
r, r« 



(1) 



=<i-fo(ri+ j-)+p(*+ -^-r-pa-fO, 

r i=V{f^f+^H\ r 2 =V(£-l+ i «) 2 +> 7 2 +£ 2 . 
The points of equilibrium are the solutions of the equations t 

There are two sets of points which satisfy these equations. One set consists 
of the two points which are at the vertices of the two equilateral triangles on 
the opposite sides of the line joining the finite bodies. The orbits which are 
asymptotic to these points and also to the periodic oscillations near these 
points are discussed in another paper. $ The other set of points consists of 

*Moulton, "Introduction to Celestial Mechanics" (1914), p. 279. 

f Moulton, "Introduction to Celestial Mechanics," p. 290; Charlier, "Die Meehanik des Himmels," 
Vol. II, pp. 102-111. 

% This paper is now under consideration for publication. 
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three points which lie on the straight line joining the finite bodies. Let the 
coordinates of these points be denoted by £„> 0, 0. where the particular value 
of £ depends upon the equilibrium point in question.* The points themselves 
will be denoted by (a), (b), and (c), where (a) lies between -fee and the 
finite mass [i, (b) between pi and 1— (i, and (c) between! — ft and — qo . 

If the infinitesimal body is given a small displacement from an equilibrium 
point and a small velocity with respect to the finite masses such that 

£ =%o+%, n =0+#, f=0 + z, 

Z>=Q+W, yi'=0+W, £=0+£', 

then the differential equations (1) become t 



(2) 



7X 

s"=-^=-=zZ(x,y 2 ,z i ), 



U =i(l^)(4+^)+^(rl+^j-^-i(i(l-(i), 



(3) 



where X, Y, and Z are power series in x, y 2 , i 2 . These series converge within 
certain regions about the equilibrium points. $ 



§ 4. The Pebiodic Obbits. 

In showing the existence of periodic solutions of the differential equations 
of the "Oscillating Satellite" which correspond to equations (3), and later in 
making the construction of these solutions, the transformations 

x—ex, y = ey, z = sz, t— 1 = (l + 5)t (4) 

are made, where e is an arbitrary parameter and $ is determined as a function 
of e so that the solutions in x, y, and z shall be periodic with the periodic 
27t in t. On denoting differentiation with respect to t by a dot over the 
variables, the differential equations (3) become as a consequence of (4)§ 

x-2(l + h)y= (l+S) 2 [X 1 +X 2e +....+X* e *-M-....], ■ 

i/+2(i+8)*=(i+a) 2 [r 1 +r 2 e+....+r jte *- i +....], \ (5) 

'£= (l+S) 2 [Z 1 +Z 2e + . . . . +Z k ^+ . . . . ], 



♦"Oscillating Satellite," equations (4). 
t " Oscillating Satellite," § 77. 



t " Oscillating Satellite," equations (6) . 
§" Oscillating Satellite," equations (11). 
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where X k , Y k , and Z k are homogeneous polynomials of degree k in x, y, and z. 
From (3) it is obvious that the X k are even in y and z, Y k odd in y and even 
in z, and Z k even in y and odd in z. The explicit values of these terms up to 
&=3 are* 

X 1 =(l + 2A)x, X 2 =iB(—2x z +y 2 +z i ), X s = 2C(2x 3 —3xy i —3xz 2 ), 
Y 1 ^(l-A)y, Y,=3Bxy, Y& =§C{-±x*y+y*+yz*), 

Z 1 =—As, Z 2 =3Bxz, Z^Ci—lrfz + ytz+z 3 ), 

- 1 ~f t i V- 

f f JO). , 



A = 



rf 3 



B 



C 



± rf> 4 



r (0)« > 
— '2 






1-f 

f <0)< 



+ 



7|oy. 



The upper, middle, or lower signs are to be taken in B according as the 
equilibrium point is (a), (&), or (c) respectively. 

If periodic solutions of (5) exist, their periods are determined from the 
periods of the solutions of the linear terms of (5). The solutions of these 
linear terms are f 

x = K x e iar + K 2 e- iaT + K B e" T + K^-" 7 , i = V^-L, 

y=in{K 1 e iaT —K 2 e- iST )+m{K s e' ,T —K i e-'' T ), 



z =K 5 cos VA'v+Ke sin VAt, n= 



g 2 +l + 2^ 
2<j ~ ; 



o 2 —l — 2A 
m = 2p 



(6) 



where K ly . . . . , K 6 are the constants of integration, and a 2 and p 2 are the 
negative and positive roots respectively of the quadratic in X 2 , 

X 4 +(2— A)tf+(1— A)(1 + 2A)=0. (7) 

There are three real periods in these solutions, viz., 2n/VA, 2it/a, and 
P=2jn/VA=2kn/a, where j and k are positive integers and VA and a are 
commensurable.^ These are the periods of the orbits of Class A, Class B, 
and Class C respectively. 

Orbits of Class A. 
The periodic solutions of Class A are § 
cc 1 =ex 1 =0e+ (ai + feicos 2s/An:)e 2 + 0s s + 

y 1 = sy 1 =0e+ (c x sin 2 VIt)£ 2 +0£ 3 + 

1 



Zi=SZ 1 : 



VI 



sin VAx ) e+0e 2 +^(3 sin V^t— sin3 V^T)e 3 + • 



$ = S 1= 0e+Sf e 2 + .... + S?'V>', 



(8) 



*" Oscillating Satellite," equations (53), (15), (35). 
t " Oscillating Satellite," equations (16), (27), and (28). 

$ In § 83 of the " Oscillating Satellite " it is shown that there are infinitely many values of p 
between and i such that ~Va and a are commensurable. 
§ " Oscillating Satellite," equations (71) . 
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where 



-35 



ttl ~ 44(1 + 21)' 



6i= 



c,= 



#»= 



—SB 



3B(1 + 3A) 
44(1— 7A + 18A 2 )' 

3 T 3B 2 (1+3A) 



V A (1-7 A + 18 A 2 )' W1 644 5 / 2 
9 T 3B 2 (1— 3A + UA 2 ) 



3 f 3B 2 (l+3A) 1 

ai ~64:AWll-7A + 18A* C J' 



[ 



16A 2 L(l-F 2A) (1—7A+18A 2 ) 



-a]. 



The expressions for x^ and ^ are power series in s 2 with sums of cosines 
and sines respectively of even multiples of ~VAx in the coefficients. The highest 
multiple of V-4t which appears in the coefficient of e k in each series is k. The 
solution for z x is a power series in odd powers of s having sums of sines of 
odd multiples of y/Ai in the coefficients. The highest multiple of VAt which 
occurs in the coefficient of e k is likewise k. The initial conditions are so chosen 
in the construction of these solutions that 

s(0)=o, i(0)=s, 

that is, the parameter e is proportional to the initial projection of the infini- 
tesimal body from the plane of motion. The actual projection is e/(l+8i). 

Numerical examples of these orbits have been considered in the " Oscillating 
Satellite " and the following results have been obtained when the ratio of the 
finite masses is ten to one, or 1 — (i=10/ll, (1=1/11 * 

TABLE 1. 



Coefficient. 


Point (a). 


Point (6). 


Point (o). 


A 


2.548 


6.510 


1.082 


a 2 


2.811 


6.820 


1.144 


P 2 


3.359 


11.330 


0.226 


ii 


2.657 


3.990 


2.014 


m 


—0.747 


— 0.397 


—3.091 


B 


6.548 


—10.961 


—1.136 


G 


18.283 


55.740 


1.196 


«! 


—0.316 


0.090 


0.249 


h 


0.151 


— 0.036 


—0.230 


Ci 


—0.112 


0.018 


0.226 


d t 


—0.037 


— 0.020 


—0.002 


$P> 


0.184 


0.467 


0.001 



* This was the ratio used by Darwin in the memoir cited at the beginning of this article. The same 
ratio was adopted by Moulton in the "Oscillating Satellites" so that he might compare his results with 
those of Darwin. 



Equilibrium Points in the Problem of Three Bodies. 



85 



(9) 



Orbits of Class B. 

The periodic solutions of Class B are * 

% 2 = e x 2 =(co& <7T)e+ (a 2 -\-b 2 cos <tt + c 2 cos 2<7T)e 2 + . . . ., 
y 2 =ey 2 = ( — wsino"r)e+ ( — nb 2 sin ffT + ^ 2 sin 2o"c)c 2 + .... 
8 =3 2 =0 <? +5 2 ( V+4 (3 V'+ , 

where a 2 , b 2 , c 2 , d 2 , hf*, . . . . , are known constants. The algebraic expressions 
for these constants may be found in equations (101), (103), and (107) of the 
" Oscillating Satellite." Their particular values for 1— |K=10/11, [i=l/ll are 
to be found in Table 2. The first seven coefficients which appear in Table 1 
have been used in determining the values in Table 2. 

TABLE 2. 



Coefficient. 


Point (a). 


Point (6). 


Point (c). 


a 2 


—4.078 


8.172 


0.554 


h 


1.996 


— 4.976 


—0.038 


c* 


2.082 


— 3.196 


—0.516 


nb 2 


5.305 


—19.855 


—0.077 


d 2 


1.250 


— 1.478 


—0.276 


hf> 


3.955 


8.553 


—1.407 



The expressions for x 2 and y 2 in (9) are power series in e with sums of 
cosines and sines respectively of or in the coefficients. The highest multiple 
of ax which occurs in the coefficient of e k is Jc. The initial conditions are so 
chosen for the construction of these solutions that 

x 2 (0)=e, y,(0)=0. 

In this case the parameter e denotes the initial displacement along the £-axis 
of the infinitesimal body from an equilibrium point. 

§ 5. The Equations op Vabiation. 

Suppose now that the initial conditions for the periodic orbits are 
changed slightly. Then the infinitesimal body will deviate from its periodic 
orbit, and the amount and character of the deviation will depend upon the 
initial conditions chosen. The equations for the disturbed orbit, in terms of 
the variables x, y, and z, will be 

a=^+p } ., y=y i +q i , z-Zj+rft, (10) 

♦"Oscillating Satellite," equations (97), (105). 

fThe r t and r % in this equation and in subsequent equations are different from r, and r 2 in 
equations (1) and (3). 
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where ex iy ey jf and eZj represent the periodic orbits of Class A or Class B 
according as j=l or 2 respectively. The additional terms p s , q h and r i are 
functions of t which depend upon the initial disturbance from the periodic 
orbit. We desire to determine these functions so that they shall have the 
asymptotic form already referred to in § 2. 

On substituting (10) in (5), the differential equations which define p f , q jf 
and r 4 are found to be 



j 



where 



p,_2(l + 8,)2,+ (l + 8,)»(P A p,+P„ 3 ,+P„r,) = (1 + 8,)'P„ 

ij,+2(i+a,)p,+ (l+^ViQflPi+Q^+Q^) = (1+ W-, \ (ll) 
r,+ (l+^nR^+R^+R^) = (l+8,)*fl„ 

P il =—l—2A+§Bex i —GCe t {2x)—y)—z)) + ...., 
P tt =Q fi =—3Bey,+12Ctx,y i + ...., 
P. 3 =P A = —ZBbz, + \2Cixfa + ...., 
Q fi =-l+A-3Bax f +lC* (4z 2 -3*/ 2 -3 2 ) + ...., 

Q iS =R i2 =—3Ce 2 y i z i + , 

R js =A— 3Bex j +kCe*(4:X i i — y)— 3s 2 ) + ...., 
^=fP e (-2^ + ^+r|) + ()e 2 +...., 
(?,=3P e p A .+ ()a 2 +...., 
P,=3P^r y +()e 2 + 

As we propose later to construct numerical solutions which are asymptotic 
to the particular periodic orbits mentioned above, we insert here the particular 
values of the various P jk , Q jh , and R jk which are obtained by making use of the 
numerical values listed in Tables 1 and 2. The following results are found. 

Orbits of Class A. 
Equilibrium point (a). 

P u =:~6.096 + e 2 (8.384— 15.686 cos 2VIt) + ()e 4 + ... ., 

P 12 =Q n =s 2 (2. 200 sin 2VZt) + ()e 4 + . . . ., 

P U = B U =— e(12.297sin V2t) + ()e s + • • • •, 

g i2 = 1.548 + e 2 (0. 833 + 2. 409 cos 2VIt) + ()e 4 +...., 

Q u =R a =0#+.... 

P 1S =2.548— ^(9. 918— 13.160 cos 2 V^t)+ ()e 4 + 
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Equilibrium point (&). 

P u = — 14. 020 +e 2 (19. 777— 23.328 cos 2VIt) + ()e 4 + ... ., 

P 12 =# 21 =e 2 (0. 592 sin 2 VlT) + (y +...., 

P 13 =R u =e (12. 890 sin VIr) + ()e 3 +. . . ., 

# 12 =5.510 + e 2 ( — 3.464 + 5.240cos2VJt) + ()£ 4 +...., 

Q 13 =P 12 =() e 8 +...., 

i? 18 =6.510+e 2 (— 16. 312+18. 108 cos 2V^t) + ()«*+ 

Equilibrium point (c). 

P u =— 3.164+e 2 (1.616— 1.746cos2V2t) + ()s 4 +...., 

P 12 =# 21 = 6 2 (0.770 sin 2VZt) + () f *+ . . . ., 

P 13 =P u=e (3.275 sin V2t) + ()e 3 + . . . ., 

# 12 =0.082+6 2 (0. 020+0. 044 cos 2VZt) + ()e 4 + . . . ., 

Q 13 =R 12 =()e*+...., 

R 13 =l. 082+e 2 (0. 020+0. 044 cos 2VAx) + 

Orbits of Class B. 
Equilibrium point (a). 

P 21 =— 6.10 + e(39.30cos<jT) + , P 22 =# 21 =<-(52.27 sin or) + , 

# 22 =1.55— e(19.64 cos an) + , P z3 =Q 23 =R n =R 22 =0, 

B 23 =2.55— f(19.64cos<JT) + 

Equilibrium point (b). 

P 21 =— 14.02— e(65.77cos<nr) + , P 22 =Q 21 = — f (131.20sin<xT) + , 

<2 22 =5.51+e(32.88cos<TT) + ...., P 23 =(? 23 =P 21 =P 22 =0, 
i? 23 =6.51+e(32.88cos(JT) + 

Equilibrium point (c). 

P 21 =— 3.16— e(6.82 cos err) + , P 22 = # 21 = 6 (6.86 sin an) + , 

Q 22 =0.08 + e(3. 41 cos ctt) + ...., P 23 =Q 23 =R a =R 2z =0 ) 

P 2S =1.08 + e(3.41cos<nr) + 

If we consider only the terms of (11) which are linear in p it q jf and r } , 
we obtain the equations of variation* For the orbits of Class A they are 

&_2(i+a l )g 1 + (i+SiWuPi+Puai+JVi) =0, " 

q 1 +2(l+$ 1 )p 1 + (l+B i y(Q u p 1 +Q 12 q 1 +Q 13 r 1 )=0, \ (12A) 

r i+ (l+^i)'(BuPi+Buqi+Buri) =0, . 
and for orbits of Class B, they are 

F 2 -2(l+S 2 )g 2 +(l+$ 2 ) 2 (P 21 p 2 +P 22 2 2 )=0, 

q\+2(l+S 2 )p 2 +(l + 8 2 y(Q 21 p 2 +Q 22 q 2 )=0, [ (12B) 

r 2 +(l+5 2 ) 2 P 23 r 2 =0. 

* Poincare, loo. cit., Vol. I, Chap. 4. 
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The periodic solutions (8) and (9) are called the generating solutions* of 
(12A) and (12B) respectively. 

§6. The Solutions of the Equations of Vakiation (12A). 

The equations of variation are linear differential equations having periodic 
coefficients, the periods being the same as those of the corresponding generating 
solutions. Differential equations of this type were first treated by Hillf in 
his celebrated memoir on the lunar theory. A very extensive list of references 
to the literature of these differential equations was given by Baker in 
a recent memoir, t 

The method which we shall adopt in the construction of the solutions of 
these equations is the one developed by Moulton and Macmillan.§ The form|| 
of the solution is, in general 

e"*(<r), 

where $(t) is a periodic function having the same period as the coefficients 
of the differential equations, and a is a constant called the character- 
istic exponent. 

It has been shown by Poincare fl that if the generating solutions contain 
an arbitrary constant which does not appear explicitly in the original 
differential equations, then a set of solutions of the equations of variation is 
obtained by taking the first partial derivatives of the generating solutions with 
respect to this constant. The characteristic exponent associated with this set 
of solutions is zero, and the solutions themselves are either periodic or consist 
of periodic functions plus i times other periodic functions. The initial 
time t is one such arbitrary constant which is present in all differential 
equations of dynamics in which the potential function e.g., U in equations (1), 
does not contain the time explicitly. Differentiation of the generating solu- 
tions with respect to this constant yields a set of periodic solutions. Another 
such constant usually present is the scale factor of the generating solutions. 
Differentiation of the generating solutions with respect to this constant yields 
a set of solutions comprised of a periodic function plus K'z, K a constant, 
times the periodic function obtained by differentiating with respect to t . 

* Poincare, loo. oit. Vol. I, Chap. 4. 

t G. W. Hill, " Collected Works," Vol. I, p. 243; Acta Mathematica, Vol. VIII, pp. 1-36. 

t H. F. Baker, " On Certain Linear Differential Equations of Astronomical Interest," Philosophical 
Transactions of the Royal Society of London, Series A (1916), Vol. 216, pp. 129-186. 

§ Moulton and Macmillan, "On the Solutions of Certain Types of Linear Differential Equations 
with Periodic Coefficients," Amebican Journal of Mathematics, Vol. XXXIII (1911), No. 1, pp. 63-97. 

|| Floquet, Annates Bcientifiques de I'Ecole Superieure, Series II, Vol. XII (1883) , p. 47. 

fl Poincare, loc cit., Vol. I Chap. 4. 
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Let us consider first equations (12A), and those solutions which have 
characteristic exponents different from zero. These equations are simul- 
taneous and must be considered together. To find their solutions let 

p 1 =JT 1 e aiT « 1 , Qi=Lie aiT v 1 , r 1 =M 1 e aiT w 1 , (13) 

where K x , L lf and M x are arbitrary constants. An existence proof, which is 
omitted here, would show that 

a x = a<°> + afV+ .... + a f ft) f 2ft + . . . ., 

u x =<>+<V+ .... +<*V*+ . . . ., 

v x = v <<»+vPe*+ .... +<*> e 2 *+ . . . . , 

w 1 =wPe+wi 3 '>e 3 + .... +wi 2k+1) s u+1 + .... 
where the af fc) are constants so determined that the various ui m , vf k \ and 
wf h+l) shall be periodic with the period 27t/VA in t. On substituting (13) in 
(12A) and considering only the terms independent of s, it is found that a{ 0) 
must satisfy the biquadratic (7) if the arbitrary constants K lf L lt and M x are 
to have values other than the trivial ones K 1 =L 1 =M 1 =0. It is shown in the 
" Oscillating Satellite," § 82, that this biquadratic admits two purely imaginary 
solutions, denoted by +<n and — <si, and two real solutions, denoted by +p 
and — p, for all values of (i such that <[i<i, and for each equilibrium point 
(a), (&), and (c). 

Since the functions u lt v lf and w x are multiplied by arbitrary constants, 
we may assume without loss of generality, that u x (0) =1. Hence 

<>(0)=1, i»P»>(0)=0 (*=1, 2, ....oo). (15) 

Then ^ x (0) and w x (0) can be determined from the differential equations which 
define v x and w x . 

The various solutions of (12A) may be constructed by substituting (14) 
and (13) in (12A), and giving to a[ 0) the values <si, — ai, p, and — p in turn. 
By equating the coefficients of the various powers of s we obtain sets of differ- 
ential equations which can be integrated step by step. The arbitrary constants 
of integration arising at each step can be uniquely determined from the condi- 
tions (15). The various a[ m can likewise be uniquely determined by the 
condition that uf h> , vf h) , and wf h+v> shall be periodic with the period 2n/'\/A 
in t. From Poincare's extension to Cauchy's theorem,* it is known that these 
solutions will converge for e sufficiently small numerically and for all values 
of t such that < i < T, where T is an arbitrary period chosen in advance. 
Obviously, we may choose the period 2n/VA. 

* PoincarS, loo. eit., Vol. I, p. 58. 
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Let us denote the solutions of (12 A) when ai 0) =ai by 

p 1 =e < "- T u u , qi=e ifflT iv n , r 1 =e i " r iw u . (16) 

The arbitrary constants are omitted here but they appear in equations (25). 
It is found that a x is a power series in e 2 with real constant coefficients. The 
functions u u and v n are likewise power series in e 2 , but their coefficients are 
sums of cosines and i times sines of even multiples of VAt. The function w u 
is a power series in odd powers of e, vanishing with e, and the coefficients are 
sums of cosines and i times sines of odd multiples of V'Az. In each of these 
three functions, the highest multiple of VAi which appears in the coefficient 
of s k is k. So far as the computation has been carried out, we have 

«ii=l+()£ 2 +.--., v u =w+()e 2 +...., 



'i 



VT 

2 r 9ff> 2 (l-13^)-(3 + 7;l-22^ 2 U W(l—A)jl+2Ay \ t 
~ <T+£ L 164a* (44— o») ? J + 



Since the differential equations ^12A) do not contain i, they are unaltered 
by a change in the sign of i. If then we change the sign of i in the differential 
equations and in the initial conditions (15), we obtain solutions which differ 
from (16) only in the sign of i. Let this set of solutions be denoted by 

Pi=e~ iffT u 12 , qi=—e- iffT iv lu r x =— e~ i<rT w 12 , (17) 

where u^ , v 12 , and w n differ from u n , v n , and w n respectively only in the 
sign of i. The same result would be obtained by putting a^ = — ia and pro- 
ceeding as in the construction of the solution (16). 

Likewise by putting <x,f> equal, first to p, and then to — p, we obtain the 
respective solutions 

Pi=e Ar « M , 2i=e PlT i>is, r 1 =e" lT w 18 » (18) 

and 

ft=e-* T «u, qi=—e~ hT v u , r 1 =—e-"' T w li , (19) 

where p x is a power series in e 2 with real constant coefficients. Further u^ , u u ; 
i'ls > v u ; w ls , w u are similar to u u , v n , and w n respectively, except that the 
coefficients of the former functions are all real ; also 

« u (t)=w 13 (— t), v u ^)=v n {~ t), w u (t)=w 13 (— t). 

The computation for the first terms of these series gives 

"is =1+ ()e 2 + , v 13 =m+ ()e 2 + , 

32? r 2 1 _ "J 

Wn ~ ~~ 44 + p 2 |_T cos VAt— -7= sin VAt J e+ ( ) /+ , 

91? 2 jp 2 (l-134)+3 + 74-224 ! f S< 2 >(l-4)(l+24)] 



Pi =P + e 2 [ 



14.4p 3 (44+p 2 ) + p 8 J + 
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Obviously, the solutions (19) differ from (18) only in the sign of t, q lt 
and r x . This property can be derived directly from the differential equations, 
which we proceed to show. We shall show that if 

Pi^FflW, ffi=^(t), r,=Ft(<v) (20) 

are solutions of the differential equations (11), then 

Pi=Ffli—t), %=—&»{—*), r,=— *•„(— t) 
are likewise solutions. If this property holds for equations (11) it will 
evidently hold for equations (12A) or (12B), i. e., when the right members of 
(11) are omitted. 

Let us consider the differential equations (1). Since U is even in v\ and £, 

it follows that -^p- is even in vj and £, -=— is odd in yj and even in £, and -~~ 

is even in v\ and odd in £. After the substitutions (2) are made, the right 
members of (3) are functions of y and z as indicated. As the substitutions (4) 
do not alter the parity of y and z, then it follows that the right members of 
(5) have the same parity in y and z that the corresponding right members 
have in y and 5. When the substitutions (10) are made in (5) we obtain 
(11), and the terms of (11) which carry (1+Sj) 2 as a factor in each equation 
arise from the corresponding right members of (5). Let these terms be 
denoted by (11, 1), (11, 2), and (11, 3) respectively. Then they have the 
following properties : 

(11, 1) is even in \y j} q^, and even in \z jf r,[, 
(11, 2) is odd in \y t , qj\, and even in \z f , r,j, 
(11, 3) is even in \y } , q s \, and odd in \z jf r t \, 

where the braces \ \ denote that the variables within are to be considered 
together. Further, since y is a factor of the right member of the second 
equation in (3), it follows that (11, 2) carries a factor which is odd in \y } , q } \. 
Similarly, since z is a factor of the right members of the third equation in 
(3), then (11, 3) carries a factor which is odd in \z t , q^. 

Now, if the signs of y n z jt q if and r t be changed in the above expressions, 
then (11, 1) remains unchanged while (11, 2) and (11, 3) change signs. On 
examining the generating solutions (8) and (9) we observe that changing the 
signs of y i and z f is equivalent to changing the sign of 1 in these solutions. 
Thus, if we change the signs of t, q } , and r,-, (11, 1) remains unchanged while 
(11, 2) and (11, 3) change signs. If these changes are made in all the terms 
of (11), the first equation remains unchanged while each term of the remaining 
equations changes sign and the factor — 1 can be divided out. Thus the 
equations (11) remain unchanged if we change the signs of t, q f and r } -. 
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Suppose now that (20) is a set of solutions of (11). Let 

Pi=P» Qt=—Qi> 7 i=— r i> ^=— T - ( 21 > 

Then the differential equations obtained by making these substitutions in (11) 
are identically the same in $ i ,q. i ,f i , and t as (11) are in the former variables. 
With the same initial conditions for p f , q } , r n and t as for equations (11), we 

obtain the solutions 

p,=F A (f), ft=F Ji (T), #v=F„(t), 

where F ilt F fi , and F jZ are the same functions as in (20). On restoring the 
former variables by (21) we have 

^=2<V(-t), q^-F^-x), r } =-F (-i) 

as solutions of (11). Therefore if a set of solutions of (11) is known, another 
set can be obtained by changing the signs of t, q } , and r f in the former set. 

Let us next consider the two remaining sets of solutions of (12A). These 
solutions have characteristic exponents zero and may be obtained, as already 
referred to in this section, by differentiating the generating solutions (8) with 
respect to the arbitrary constants t and s. Taking the constant t first, we 
have the solutions 



*=f = § £~T^[<**V2«.IV2«>*.<>*+. 



'■=^4l=-mh v3 " t+ w+->' 
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1+Sx 
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1+Si 


1 






-*»»• 



%5, 



As these solutions are later multiplied by arbitrary constants, we may neglect 
the factor — z , . and consider 

1 + 01 

Pl = «15f fll = «16» *1 = W« (22) 

as the solutions. Similarly, on differentiating (8) with respect to e we obtain 

dxi dy! dsx 

Pl= a7' qi= -df> ri= ~d7- 

Since t= ($—*„) /(l+o\) and ^i is a function of e, then the constant e enters 



A 1= 
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the generating solutions not only explicitly but implicitly through 5i and t. 

Hence 

dxx _ (SxA , S%i dx_ d$i Sy^ _ (SyA 3|/i <fo_ d$± 
Se ~ \ Se ) + 3t Si,. Se ' Si ~ \ Se I + 3x 35, Se ' 

35, _ /35A , 3«i 3t 35j 
"37 _ V5T7 + ~3T W, "37' 

where the parentheses around the partial derivatives denote that the differ- 
entiation is performed only in so far as e occurs explicitly. When the 
differentiations are performed, we have 

Pi=v u +E^u 15 , q!=v lf) +Ki:v 1B , r 1 =w M +Kr:w ls , (23) 

where 

«i.= (-j| 1 ) = 2(a 1 +& 1 cos2V2T)£+()e 8 +...., 
Vw = (l?)= 2(Cl 8in2 VlT)«+ ()e 3 + . . . ., 

"i"/" v2 sin VjT+()f2+ * • " '' 

^ = - T ^-[23? , E+() e 3 +....]. 

It remains now to show that the solutions which we have obtained, viz. 
(16), (17), (18), (19), (22), and (23) constitute a fundamental set. The 
criterion for a fundamental set is that the determinant formed from these 
solutions and their first derivatives with respect to t shall be different from 
zero. This determinant is 
e" ,T M u , e- to,T u u , e'' r « u , e~* ,T WM> w 15 , u^+Knu^, 

e i " T (w 1 M 11 -(-»< 11 ) ) e-*" T (— i<Ti«i 2 +Mi 2 ), e" r (p 1 M 1 3+M 13 ), e-" T (— pi«u+«i 4 ), % 5 , «M+JC(tM«+«u)> 
e""w u , e-*" T (— iv n ), e hr v ia , e-" r {—v u ), v w ,v w +K'iv li , 

«""(-%+»%), e-*" T {c 1 Vu—iv vl ), e^^Vu+bu), e-* iT (p 1 v u — v u ), i> u , v u +K(ti) u +v u ), 
e" lT iw a , e-*' iT (—iw n ), e"' T w a , e- hT (—w u ), w u , w 1<s +Ki:w 16 , 

e' VlT (— SiW n +iw n ), e-"' T (<7 1 w n —iw lt ), e'^^w^+w^), e-' lT (p 1 Wi 4 — w u ), w^^+K^w^+w^). 

It is a constant,* and its value can be determined with the least difficulty by 
putting x=0. Thus we obtain 

A 1 =iis(mp+na) (ma— no) + terms of higher degree in e. 

It is shown in the "Oscillating Satellite," equations (36), that mp+na and 
ma — np are different from zero. Hence the determinant Aj is different from 

♦ Moulton, " Periodic Orbits," Chap. 1, Sec. 18. 



(24) 



(25) 
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zero for e not zero, but sufficiently small numerically. Therefore the solutions 
which have been obtained for (12A) constitute a fundamental set, and the 
most general solutions are 

p 1 =K 1 e ic * T u n +K 2 e-^ T u n +K z e"- T u li +K i e~ hr u u +K & u is 

+K 6 (u 16 +Ktu Vi ), 

q 1 =K 1 e i ' lT iv 11 —K^-^iv^ +K 3 e hT v 13 — Z 4 e-"' T v u +K 5 v 15 

+K 6 (v 16 +K'zv 15 ), 

r 1 =K 1 e i ' lf iw n —K 2 e- i '' T iw n +K a e'' lT w 13 —K i e-'' lT w u +K 5 w 15 

where K lf . . . ., K 6 are arbitrary Constants. 

§7. The Solutions of the Equations of Vabiation (12B). 

The solutions of the equations (12B) are obtained in the same way as the 
preceding solutions were found. The construction is simplified by the fact 
that the last equation is independent of the first two. 

The two sets of solutions of the first two equations of (12B) which have 
characteristic exponents different from zero, are 

p 2 =e"' r u n , q 2 =e>* r v 21 , (26) 

and 

p 2 =e-"' r u 2i , q 2 =—e-"* T v 22 , (27) 

where o 2 is a power series in e with real coefficients, the terms of lower degrees 

being 

p 2 =p+0e+p|V+p£V+ 

The functions w 21 , v 21 , u 22 , and v 22 are power series in e with sums of sines and 

Cosines of multiples of at in the coefficients, the highest multiple of at which 

occurs in the coefficient of e* being k. All the numerical coefficients in these 

series are real. From the property of the differential equations (11) which 

was proved in the preceding section, it follows that the solutions (27) are 

obtained from (26) by changing the signs of t and q 2 . Thus w 22 and v 22 differ 

from m 2 i and v 2X respectively only in the signs of the sines. 

Since these solutions (26) and (27) are multiplied later by arbitrary 

constants, we may take w 21 (0)=l. The algebraic forms of the functions u 21 

and v a are 

*%= 1 + (a 21 cos at + b 21 sin <jT)e-f ()e 2 + , 

v 21 =m+ (o 22 cos <7T+ 6 22 sin<7T)e-j-()e 2 +. . . ., 

m 22 (t)=m 21 (— t), v 22 (t)=v 21 (— t). 
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The two remaining solutions of the first two equations of (12B) are 
obtained by differentiating the generating solutions (9) with respect to t and s. 
Differentiation with respect to t gives 

9% _ 1 dx 2 _ dy 2 _ 1 dy 2 
Pt ~ dt ~ l+3 2 8t ' 22 ~ dt ~ l+3 2 3t ' 

and, since the multiplier — t may be dropped, the solutions become 

l + 2 



p 2 =u 23 = -^- =( — a sin<TT)e — (b 2 a sino , T + 2c 2 <7 sin 2<tt)c 2 + . . . 
q 2 =v^ = -Jr- = ( — na cos <tt)6 — (nb 2 cr cos <rv — 2d 2 a cos 2<TT)e 2 + 

<7T 



(28) 



Differentiation with respect to e gives 

P2=u 2i +L%u ii , q 2 =v 2t +LTV 23 , (29) 

where 

u 2t = (-J^)= cos<TT + 2(a2+& 2 cos<rir + c 2 cos 2<rir)6+ , 

v 2l = (-Jp) = — n sin crv— 2(«& 2 sin ctt— i 2 sin 2<nr)e+ , 

L = -TXF (2S| 2 > e +3SfV+. .. .)• 
1 + &2 

The determinant of these solutions (26), (27), (28), and (29) and their 
first derivatives with respect to t is 

e p ' T u n , e- p ' T u n , «28, m 24 +Ltm 28 , 

e" 3T (p2U a +u 21 ), e- p ' T (— p 2 u n +u 22 ), i%, u^+L^u^+u^), 
e fi ' T v 2l , e~' sT (— v 22 ), v u , v u +Lwu, 

e p ' T (p 2 v 21 +v 21 ), e~ HT {p 2 v 22 —v n ) , t> 28 , Vu+L(tVa + v u ), 
= 2<js(mp+n<s) {ma— rip) + terms of higher degree in e. 

As in the former determinant A 1} equation (24), the factors mp+na and 
m(T _ M p are different from zero, and since <r=£0 it follows that A^O for e not 
zero, but sufficiently small numerically. Hence the solutions (26), (27), (28), 
and (29) constitute a fundamental set of solutions of the first two equations, 
and the most general solutions are 

p 2 =L 1 e^ T u 21 +L 2 e-''' T u 22 +L s u 23 +L i (u 2i +LT!u 2S ), "I 
q 2 =L 1 e"' T v 21 — L 2 e- psT v^+L 3 v 23 +L i (v 2i +L'vv 23 ), J 

where L x , . . . . , L 4 are arbitrary constants. 
13 



A,= 



(30) 
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The general solution of the last equation of (12B) is readily obtained, 

and it has the form 

r 2 =L 5 e iaT w Z i+L 6 e- iaT w 22 , (32) 

where L 6 and L 6 are arbitrary constants, and o, w 21 , and w 22 are power series 

in s of the form 

3Be — 

w 21 =l-\ — 27-5 — j-tt- (<r 2 cos <tt — 2i\ 'A sin a% — a 2 ) + ...., 

3Bs „ . 1 — 

w 22 =l-) — 2 2 j-rr- (<rcos <JT+2tVJ. sinaT— a 2 ) + ...., 



6)=V^ + 0e+« 2 e 2 + 



Since L s and L 6 are arbitrary, the initial values of w 21 and w 22 may be chosen 
so that w n (0) =w 22 (0) =1. The coefficients of the various powers of s in these 
functions are sums of cosines and i times sines of multiples of <tt, the highest 
multiple which occurs in the coefficient of e h being k. Further, w a and w 22 
differ only in the sign of i. In the series for a, all the coefficients are real. 
The determinant of the two solutions in (32) and their derivatives with 

respect to t is 

e i0T w 21 , e- iaT w 22 , 

e iaT {iaw 2l -\-w 21 ) , e~ iaT (— iaw 21 +w 21 ) , 
—2ia-\- terms in e, 



A a = 



and this is different from zero for | e | sufficiently small. Hence the solutions 
(32) constitute a fundamental set of solutions of the last equation of (12B). 

§8. The Consteuction of Asymptotic Solutions. 

Having determined the solutions of the equations of variation, we shall 
now proceed to construct the asymptotic solutions of (11). "We shall consider 
first the solutions which approach zero as 1 approaches +oo, and then show 
how to obtain from these the solutions which approach zero as t approaches 
— 00. 

In making the constructions, it is convenient to introduce a parameter y 
by the substitutions 

Pi=PiY> 2,= <Z/7> r i=W> ( 33 ) 

where p } , g,- , and r i are new dependent variables. Then, since the asymptotic 
solutions are to be constructed as power series in y, we put 

P i =PP+pPr+pfV+-..-, 2,=sf+<zfV+2fy+ 1 m , 

r.=rf>+rPy+r?Y+ J { ' 
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When (34) and (33) are substituted in (11), the factor y will divide ont. Let 
the resulting differential equations be denoted by (11'). On equating the 
coefficients of the various powers of y in these equations, we obtain sets of 
differential equations which define the functions pf, g, w , and rj k) in (34). In 
order that the solutions shall be asymptotic we must impose the condition 

(d), that each term of the solutions for the various p\ h) , q$" y , and r, (fc) shall 
contain e~ PlT as a factor. 

As an arbitrary constant arises at each step of the integration we may impose 
the further condition 

(C 2 ),thatp i (0)=a j , 

from which it follows that 

^ 0) (0)=a,., *>H0)=0, (^ = 1, ....oo). 

Asymptotic Orbits of Class A. 

Let us consider first the orbits which are asymptotic to the periodic orbits 
of Class A. We put j=l and consider the various differential equations 
obtained by equating the coefficients of the same powers of y in (11'). 

For the terms in pf>, qf\ and rf> we obtain a set of differential equations 
which are the same as (12A), except for the superscript 0. The general solu- 
tions are the same as (25), and when the condition (d) is imposed all the 
constants of integration must be zero except the one associated with e~ PiT . 
From condition (C 2 ) it follows that this constant must have the value a x . 
Therefore the desired solutions are 

pf> =a 1 e- hr u u , d 0) = -aie- p,T Vi 4 , rf =— ai e- ftT w M . (35) 

The differential equations which are obtained from equating the coefficients 
of y to the first degree in (11') have the same left members as (12A), except 
for the superscript 1 on the variables. On denoting the corresponding right 
members by Pf», Q™, and RP, we find that 

P? = ea\e- thT U®, Q?> = eaJe-^Fff , B« = ea!e-^ T TF^, 

where TJ$, V$, and W$ are power series of the same form as u u , v u , and w u 
respectively. The complementary functions of these differential equations 
are the same as (25), but let us suppose that the arbitrary multipliers are 



(36) 
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kP, , k£\ The particular integrals can be obtained by varying these 

multipliers. Thus 

Af V fflT « u + JB^er^Uu + lcPe hT u 13 + H»e-* T u u +U iy u 15 

+ kP(u M +K™ u )=Q, 
k[ 1) e i ^ T (i<j 1 u 11 + «u) — 4 (1) e _<<riT (^Mm— « 12 ) + *J»e* r ( 9l u ls + w 18 ) 

-^e-^Cp^u-Mu) + /#>% 5 + &« [«16 + #(% 6 + ™ 15 ) ] =P?>, 

i^Mvn- A^e-^Wu+J^V^s— k^e^ T v u +'k^v 15 

+*f ) (w w +JCTW M )=0, 

^^(-qwu+wa) iiife^^^tta-Hpa) +jf ) e* r ((» l w B +w») 

+k^e-"^(p 1 w M -w u ) +iP%+4 a) [%+I(%+T%) ] =E 1 a >. 

The determinant of the coefficients of k^, , ty is A lt equation (24), and 

since it is different from zero for e not zero but sufficiently small numerically 
these equations can be solved for k^, . . . . , k£ y . Thus 

Afl) 

W=^> (1=1,.... ,6), (37) 

where AJP is the determinant formed by replacing the elements of the Z-th 
column of A lf with 0, P?>, 0, Q?, 0, and Bf> respectively. Since P®, Qj», and 

B'i i:> contain no terms in e ±iff]T or e ±hT , the integrations of (37) for k?\ , k?> 

will yield no terms in t explicitly. Such terms will occur, however, in the 
integrations for k^ and k^\ but when they are substituted in the complementary 
functions they will cancel off. The complete solutions are thus found to be 

+KV(u 16 +Kzu 15 ) +ea?e- 2 *X?, 
gF=K?& r iv u —KFe-«'>- r w a +KFe* r v n —K2 i e- hr v u +KPv» 

+KP (v 16 +Kw w ) +eaie-fv®, 
ri l) =K^e ia ^w 11 —K^e- ia ^iw li +K^e hT w Vi —K2 ) e- hT w u +K^w 15 

+KP {w u +Kvw 15 ) +ea.fe- 2l » T w$, 

where K^, , K^ are the constants of integration, and u$, v$, and w$ 

are functions similar to u u , v u , and w u respectively. On imposing condition 
(Ci) upon these solutions we obtain 



(38) 
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and from condition (C 2 ) we have the remaining constant 

Jtf> = -«afcS>(0). 

When the constants of integration are thus determined, the solutions (38) 
take the form 

p™=ea\(e^u™ +e- 2 " T «#), qP = e<x*(e-'>M\ > +e-** T t>8 > ), 1 

rf > = ealie-^wft + e~ thr w^ ) , J ^ ) 

where the functions u[?, u§ ; vff, v$ ; w[\\ w$ are of the same form as u u , 
v u , and w u respectively. 

The remaining steps of the integration can be carried on in essentially 
the same way. By an induction to the general term we shall show thai the 
solutions of (11) can be constructed to any desired degree of accuracy. 

Let us assume that p[ r) f qi" y , r[* have been constructed for v=l, 2, , A;— 1, 

and that 

f + 1 v+l v+l 

2>w=eX +1 Se-' ftr « 1 ( r ) , qP=e y a y i +1 2e- lhr vtf, rf ) =«'a[ +x S e~ lhr w^, (40) 
j=i i-i i=\ 

where the functions wfl , v§>, and wff are similar to w u , v u , and w u respectively. 
We propose to show that the solutions for p[ k) , qi k \ and r[ k) are the same as 
(40) if v=k. 

Consider the set of differential equations obtained by equating the 
coefficients of y* in (11'). The left members are the same as (12A), except 
for the superscript k on p lt q lf and r t . Let the corresponding right members 
be denoted by P[ k \ Q?\ and B?\ Then it is found that 

Pf > = s k ai +1 "Se- ,hT Ui k \ Qi» = e fc aj +1 sV'^Ftf >, flj*> = e V +1 sV'^TF^, 
i=i j=i i=i 

where Ui ky , V[ k \ and W[ k) are functions of the same form as u u , v u , and w u 
respectively. The complementary functions of these equations are the same 
as at each preceding step of the integration, and the complete solutions may 
also be obtained in the same way as equations (38) were determined. The 
equations analogous to (36) and (37) differ from (36) and (37) respectively 
only in the superscript k instead of 1. The complete solutions will therefore 
have the same form as (38), except that the respective particular integrals 
are functions similar to P[ k \ Q[ k \ and JB|*> instead of P? 5 , Qjp, and JJJ as in 
(38). From condition (Ci) the arbitrary constants will all be zero, except 
K ik) , and from condition (6' 2 ) this constant is found to carry the factor 6*aJ +1 . 
Hence the desired solutions for pi k) , q{ k) , and r\ k) are the same as (40) if v—k. 
This completes the induction. 



100 



Buchanan: Asymptotic Satellites near the Straight-Line 



, co) are substi- 



(41) 



When these solutions for p[*\ gf>, and rf°, (v=l, 2, . 
tuted in (34) and (33), we obtain as the solutions of (11), 

Pi = I 2 e-*» T <> (t ) e>- W, 
g x = 2 2 e-^ T < (t ) e>- W, 

f!=H «-***«# (t ) e>- w, 

where the superscript on the functions w$?(t), w$?(t), and w$(t) has been 
made to conform with the corresponding power of y. Since the two arbitrary 
parameters a x and y occur only in products as indicated, we may suppress 
either without loss of generality. Let us suppose that 0^=1. 

Equations (41) are the asymptotic solutions of (11) which approach zero 
as i and therefore t approach +co. In the same way there could be obtained 
the asymptotic solutions which approach zero as t and therefore t approach 
— oo. It is not necessary to construct these solutions, however, since they can 
be obtained directly from (41), as we proceed to show. 

It was shown in the two paragraphs following equation (20) that if a set 
of solutions of equations (11) is known, then another set can be obtained by 
changing the signs of t, q x , and r t in the former set. On making these changes 
in (41) we thus obtain the asymptotic solutions of (11) which approach zero 
as t approaches — oo, viz., 



ft = + E 2 e** T «&> ( — t ) e*-y , 
1=1 k=i 

& = — 2 2 e* piT «ff ( - t ) e>-y , 
j=i k=i- 

fi = — S i e khT w$ ( — x ) s'-y . 

j=l k=\ 



(42) 



Upon substituting (41) and (42) in (10), and returning to the original 
variables £, yj, £ through the substitutions (4) and (2), the asymptotic solutions 
of the original differential equations (1) in terms of i are found to be 



£=£o+»i+ 2 2 e- 4 ^ «tf ( ± <c)*y, 

;=1 fc=l 

^=0 +y x ± 2 2 e**"' 7 vff ( ± T)eV, 

i=\ k=l 

f=0 + s x ± 2 2 e T ** T w$ ( ± f ) eV- 



(43) 



Where the double signs occur, the upper signs give the solutions which 
approach the periodic orbits as t approaches +oo, and the lower signs give 
the solutions which approach the same periodic orbits as t approaches — oo. 
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Asymptotic Orbits of Class B. 

We consider now the orbits which are asymptotic to the periodic orbits of 
Class B. The method of constructing these orbits is entirely similar to the 
one used in constructing the asymptotic orbits of Class A. 

On putting j=2 and considering the differential equations obtained by 
equating the coefficients of the various powers of y in (11') , that is, in (11) 
when (33) and (34) have been substituted and y divided out, we obtain solu- 
tions which will give the asymptotic orbits, provided the arbitrary constants 
of integration are so chosen that conditions {C x ) and (C 2 ) are satisfied. 

The differential equations arising from the terms independent of y in 
(11') are the same as (12B), except for the superscript on p 2 , q 2 , and r 2 . 
Their solutions are therefore the same as (31) and (32). "When the conditions 
(Ci) and (C 2 ) are imposed, we obtain 

pf> = a 2 e-* T « 22 , gS» = -a^e-^Vn , r 2 (0 > = 0. 

From the coefficients of y to the first degree in (11') we obtain the differ- 
ential equations which define p£\ qi X) , and r£\ These equations have the same 
left members as (12B), except the superscript 1. Let the respective right 
members be denoted by P| J) , Q£\ and Rg*. Then 

pv=eale- 2 >° T Ug\ QP=ea\e-*<>* T Vg\ R?>=0, 

where U$ and V$ are functions of the same form as u 22 and v 22 respectively. 
The complementary functions of the differential equations are the same as 
(31) and (32), but let the arbitrary constants be denoted by li iy , . . . ., li iy . By 
using the method of the variation of parameters to obtain the complete solu- 
tions we have 

iPe^Un+iPe-i'Uu+iPua+iP (Uu+btUu) =0, 
^P ) e* T (p»u«i+«a) +^ 1) e- p2T (— o 2 u^-\-u 22 ) +4 (I) « 2S 

+ii l >[« 24 +L(« 23 +TO 23 ) ] =p a (1) , 

i^VK-iPe-^Vn+iUvK+ip (Vu+Dtv*) =0, 
iPeffaVn+Vn) +ipe-'* T {o 2 v 2t -v 22 ) +#>%, 

ipe {ar w 21 + iPe- iaT w 22 = 0, 
il»e iar (mw 21 +w 21 ) —iPe- iar (iow n —w 22 ) =0. 

The last two equations are independent of the first four. Since the deter- 
minant of the coefficients of if? and Itf* in these last two equations is different 
from zero, viz., A 3 , the only solutions for IfP and Iff* are l£ y =l£ y =0. 



[ (44) 
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The determinant of the coefficients of Zjp, . . . . , ip is different from zero, 
viz., A 2 equation (30), and therefore the first four equations of (44) can be 
solved for jf >, , 'Vp. The resulting solutions are 

4° = ^, (*=1,....,4), (45) 

where A$ is the determinant obtained by replacing the elements of the k-th 
row of A 2 with 0, P^, 0, and Q^ respectively. Since the right members P| w 
and QP do not contain terms in e ±p * r , the integrations of (45) for 1{ X) and l^ 
will not contain terms in t explicitly. Such terms arise, however, in the 
integration for IJP and IP, but when they are substituted in the complementary 
functions they cancel off. The complete solutions of the differential equations 
at this step are thus found to be 

pP =L^e"' T u 21 +L^e-^ T u 2Z +Li 1) u 23 +L^(u u +L'cu 23 ) +«aJe-«* T «g> (t) , " 

q™ =LPe^v 21 -LPe->* r v 22 +LPv 23 +LP(v 2l +L'zv ss ) +Eate~ 2l » T vg>(%), I (46) 

rj» =Li 1 >e i ° r w n +Li»e- iaT w 22 , 

where L^ . . . ., L^p are the constants of integration, and u${i) and v$(?z) 
are functions similar to u 21 and v 21 respectively. From condition (Cj) we have 

LP=L^=L«=L«=L^=0, 

and from condition (C 2 ) 

Lf) = ~«aSi4i'(0). 

Then the solutions (46) become 

pP = eaf [e-»M? (t ) + «-*» T iifi> (t) ] , 
qP = eai [ e~* T t;ff (t ) + e~^v^ (t ) ] , 
r 2 « =0, 

where ^(t) and #$ (t), (&=1, 2), are similar to w 21 and v 21 respectively. 

The remaining steps of the integration can be carried on in the same way. 
Proceeding by induction to the general term, we find that 

p$-> = s v a v 2 +1 2 e-^ T M$ (t ) , q™ = s"(4 +1 2 e-^vg (* ) , *f> = 0, 
fc=i fc=i 

where the functions m|J?(t) and v|?(t) have the same form as w 21 and t> 21 

respectively. When these terms are substituted in (34) and (33), we find the 

solutions of (11) to be 

p 2 =2 ile-'N^^^'al/, ft =2 Se-*^E>(T) e ^ 1 a 2 y'", r 2 =0. (47) 

As in the solutions (41) we may put a 2 =l. These equations (47) are the 
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asymptotic solutions which approach zero as t approaches +<»• By changing 
the signs of i and q % as in equations (42), we obtain 

p 2 =+2 il6** T «a>(— T)e^-y, q 2 = — 2 Se^'f-^^y, r 2 =0, (48) 

which are the asymptotic solutions approaching zero as t approaches — oo. 
Since r 2 =0 in (47) and (48), the asymptotic orbits are of two dimensions and 
are complanar with the periodic orbits of Class B. 

In terms of the original variables £, vj, £, these asymptotic solutions of 
equations (1) are 



-oo / _ 

>7=o +# 2 ± s s e+** T va> ( ±T)«y, 

7=1 ft=l 



(49) 



These solutions approach the periodic orbits of Class B as t approaches + oo 
or — oo according respectively as the upper or lower signs are taken. 

§ 9. Geometrical Considerations. 

The asymptotic solutions (43) and (49) contain the three undetermined 
constants t , s, and y. The constant t denotes the initial time and may be put 
equal to zero without loss of generality. The parameters s and y are the 
respective scale factors of the periodic and asymptotic orbits. The physical 
interpretation of the parameter e has already been referred to in the latter 
part of § 4. From the way in which the initial conditions were chosen in (C 2 ) , 
it is evident that ey denotes the ^-component of the infinitesimal body's initial 
displacement from the periodic orbit. 

Let us now consider the directions in which the asymptotic orbits approach 
the periodic orbits. We shall discuss only the orbits which approach the 
periodic orbits of Class A as t approaches +co, since the discussion is 
entirely similar for the other asymptotic orbits. 

As t becomes very large, the most important terms of the solutions of 
(11) are those in e~ hT . These terms arise only from the complementary 
functions of the differential equations which define the various pi k \ q[ k \ and 
r[ k) in (34). Neglecting the explicit values of the constants of integration 
which are associated with e~ p,T at each stage of the integration, we find that 
the predominating terms of the solutions as t approaches + co are 
p 1=+e -ps Uu[K (» r+K ay +K vy + _ . >L 

q 1 = -e-'"v u [Ki»y+K?y+K?y+ ....], 
r 1 =-e-* T w u [KPy+K?y+K?y+.. . .], 

where KP, K'p, . . . . , are the constants of integration. 
14 
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The projection on the £>7-plane of the direction of the approach is 
limit d 4i/ d Pi = limit i^IlMii. 

This limit is independent of y, but, since u u and v u contain sines and cosines 
of multiples of VAz, it is indeterminate. The projections on the ££- and 
>?£~planes are likewise found to be indeterminate. Similar results are obtained 
for the orbits which approach the periodic orbits of Class A as t approaches 
— oo, and also for the corresponding orbits of Class B. 

§ 10. Illtjstbative Examples. 

We shall conclude this article with numerical examples and diagrams of 
the periodic and asymptotic orbits which have been discussed in the preceding 
sections. In these examples the ratio of the finite masses is ten to one, or 
1 — ^=10/11 and (i=l/ll f being the ratio used in the particular periodic 
orbits already mentioned in § 4. 

The solutions for the asymptotic orbits have been carried out to the third 
degree in s for orbits of Class A and to the second degree in e for orbits of 
Class B, but only to the first degree in y for both classes of orbits. 

In the numerical results that are to be found in the tables which follow, 
e has the value 0.5 for Class A and 0.01 for Class B while ^=0.1 for both 
classes. The values of eXj, sy t , and eZj {) = !, 2) in the tables are the 
coordinates for the periodic orbits for the various values of t indicated. The 
values for eyp n eyq^, and eyr s are the amounts which must be added to the 
coordinates for the periodic orbits in order to obtain the asymptotic orbits. 
For y = 2, z j =r j =0. 

The periodic orbits are represented by the heavy lines in the diagrams, 
and the asymptotic orbits by the dotted lines. The arrows indicate the 
direction of motion. No arrows appear in the periodic orbits in Figs. 2, 5, 
and 8, as in these projections the infinitesimal body oscillates up and down 
along the same curve. The origin of coordinates is taken at the equilibrium 
point marked in the diagram, and the axes are parallel to the rotating %v;£ axes 
(see equations (2)). The unit of measurement is indicated in each diagram. 

The following results have been obtained. 

Orbits of Class A. 
Equilibrium point (a). 

ep 1 = e-* T y[e + e»(— 0.29 + 0. 29 cos 2VZt — 1.08 sin 2V2t) + ....], 

eq 1 = e-* r y[— 0.75 e+e 3 (0. 22— 0.61 cos 2VZt— 0.56 sin 2VZt) + ], 

er 1 =er* T y[f (1.58 cos VZt + 0.91 sin V2t) + ], 

p x =1.83— 0.23e 2 +...., .4 = 2.548. 
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TABLE 3. 

6=0.5, y=0.1. 



105 



T 


till! 


«!/i 


est 


*7Pi 


nil 


ryi 





— .041 








.0500 


— .0420 


.0400 


.1 


— .043 


— .008 


.043 


.0385 


— .0369 


.0359 


.2 


— .047 


— .015 


.085 


.0295 


— .0309 


.0315 


.3 


— .054 


— .021 


.125 


.0223 


— .0259 


.0270 


.4 


— .064 


— .026 


.163 


.0173 


— .0212 


.0227 


.5 


— .074 


— .028 


.198 


.0136 


— .0173 


.0187 


.6 


— .084 


— .028 


.228 


.0111 


— .0142 


.0155 


.7 


— .095 


— .025 


.254 


.0093 


— .0107 


.0118 


.8 


— .104 


— .021 


.274 


.0080 


— .0083 


.0091 


.9 


— .111 


— .015 


.289 


.0071 


— .0064 


.0067 


1.0 


— .116 


— .007 


.297 


.0065 


— .0048 


.0048 


1.2 


-.115 


.009 


.296 


.0056 


— .0029 


.0020 


1.4 


— .103 


.022 


.271 


.0045 


— .0020 


.0003 


1.6 


— .083 


.028 


.223 


.0034 


— .0015 


— .0005 6 


1.8 


— .062 


.028 


.157 


.0025 


— .0013 


— .00091 


2.0 


— .046 


.014 


.078 


.0017 


— .0011 


— .0009 5 


2.4 


— .048 


— .016 


— .092 


.0005 9 


— .0006 


— .0006 5 


2.8 


— .087 


— .027 


— .233 


.0002 3 


— .0002 7 


— .0003 


3.2 


— .116 


— .006 


— .298 


.00014 


— .0000 9 


— .0000 9 


3.6 


— .102 


.022 


— .269 


.0000 92 


— .0000 7 


— .0000 4 


4.0 


— .060 


.024 


— .151 


.0000 48 


— .0000 3 


— .0000 18 


4.4 


— .041 


— .003 


.014 


.0000 15 


— .0000 2 


— .000016 



The projections of the above orbits on the coordinate planes are given in 
Figs. 1^2, and 3. 

z 



IOJ 



Fig. 1. 




Fig. 2. 




«> X 



Fig. 8. 



106 Buchanan: Asymptotic Satellites near the Straight-Line 

Equilibrium point (b). 

ep 1 =e-" T y[e+s 8 (0.222— 0.222 cos 2VZt— 0.454 sin 2VZt) + ], 

£3l =e- p ' r y[0.40s+« 2 (0.089— 0.182 cos 2 VZx— 0.037 sin 2V2t) + ], 

sr x -e-<" T y[e(— 0.345 sin VZt— 0.523 cos VZt) + ], 

Pl =3.366— 2.837 « 2 +...., .4 = 6.510. 

TABLE 4. 
e=0.5, y=0.1. 



T 


e#, 


«2/i 


es, 


«7Pi 


«79i 


rm 





.0135 








.0500 


.0188 


— .0131 


.1 


.0147 


.0022 


.0474 


.0300 


.0119 


— .0094 


.2 


.0178 


.0038 


.0917 


.0186 


.0073 8 


— .0062 6 


.3 


.0221 


.0045 


.1355 


.0118 


.00517 


— .0038 6 


.4 


.0266 


.0040 


.1602 


.0077 9 


.0034 5 


— .0022 6 


.5 


.0300 


.0025 


.1799 


.00519 


.0022 7 


— .0009 57 


.6 


.0315 


.0004 


.1878 


.0034 9 


.0014 7 


— .0005 76 


.7 


.0307 


— .0019 


.1837 


.0023 1 


.0009 36 


— .0002 24 


.8 


.0278 


— .0036 


.1677 


.0014 9 


.0005 74 


— .0000 454 


.9 


.0236 


— .0045 


.1594 


.0009 34 


.0003 47 


+ .0000 350 


1.0 


.0206 


— .0044 


.1047 


.0005 83 


.000213 


+ .0000 616 


1.2 


.0136 


— .0007 


.0150 


.0002 04 


.0000 762 


+ .0000 497 


1.4 


.0166 


.0034 


-.0784 


.0000 750 


.0000 309 


+ .0000 363 


1.6 


.0252 


.0043 


— .1519 


.0000 308 


.0000 136 


+ .0000 0936 


1.8 


.0312 


.0011 


— .1867 


.0000137 


.0000 0590 


+ .0000 0257 


2.0 


.0289 


— .0032 


— .1837 


.0000 0594 


.0000 0232 


+ .0000 0056 6 


2.2 


.0205 


— .0044 


— .1169 


.0000 0233 


.0000 00810 


— .0000 00216 


2.4 


.0140 


— .0014 


— .0301 


.0000 0086 6 


.0000 0030 8 


— .0000 0002 2 


2.6 


.0156 


.0029 


.0645 


.0000 0029 9 


.0000 0012 4 


— .0000 00019 



The projections of the above orbits on the coordinate planes are given in 
Figs. 4, 5, and 6. 

z 



<b 



<h> 







Fig. 5. 



Fia. 4. 

Equilibrium point (c). 

ep 1 = e-'* T y[e+e s (— 0.104+0.104 cos 2 VZt— 0.817 sin 2 VZt) + ], 

e 2l -e' hT y [3 .08 e+i?(— .320— . 787 cos 2 VZt— . 106 sin 2 VZt) + ] , 

eri - e -^ r y[e(— 0.119 sin VZt— 3.145 cos VZt) + ], 

Pl =0.457— 0.025 6 2 +...., ^ = 1.082. 
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TABLE 5. 
£ =0.5, y=0.1. 
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T 


e&i 


«2h 


e«i 


£7Pi 


«79i 


tyn 





.0050 








.0500 


.1402 


— .0786 


.2 


.0099 


.0228 


.0989 


.0418 


.1279 


— .0733 


.4 


.0238 


.0418 


.1939 


.0349 


.1180 


— .0657 


.6 


.0443 


.0536 


.2803 


.0297 


.1098 


— .0560 


.8 


.0678 


.0563 


.3547 


.0264 


.1030 


— .0455 


1.0 


.0905 


.0493 


.4138 


.0245 


.0960 


— .0346 


1.2 


.1084 


.0340 


.4550 


.0236 


.0893 


—.0239 


1.4 


.1185 


.0129 


.4766 


.0234 


.0825 


— .0139 


1.6 


.1190 


— .0105 


.4781 


.0233 


.0754 


— .0051 


1.8 


.1099 


— .0320 


.4584 


.0229 


.0680 


.0026 


2.0 


.0927 


— .0481 


.4190 


.0221 


.0612 


.0088 


2.2 


.0703 


— .0560 


.3619 


.0209 


.0543 


.0135 


2.4 


.0466 


— .0543 


.2890 


.0191 


.0481 


.0168 


2.6 


.0256 


— .0434 


.2035 


.0169 


.0426 


.0187 


2.8 


.0110 


— .0251 


.1094 


.0147 


.0377 


.0194 


3.0 


.0051 


— .0025 


.0105 


.0123 


.0341 


.0191 


3.4 


.0220 


.0400 


— . 1843 


.0086 


.0289 


.0161 


3.8 


.0653 


.0564 


— .3475 


.0064 5 


.0251 


.0113 


4.2 


.1069 


.0359 


— .4517 


.0057 3 


.0218 


.0060 


4.6 


.1194 


— .0080 


— .4790 


.0056 2 


.0184 


.0014 


5 


.1101 


— .0318 


— .4723 


.00510 


.0151 


.0015 


5.4 


.0491 


— .0549 


— .2971 


.0046 6 


.0119 


— .0040 


5.8 


.0121 


— .0272 


— .1195 


.0035 9 


.0093 


— .0047 


6.2 


.0081 


.0182 


.0787 


.0014 5 


.0076 


— .0044 



The projections of the above orbits on the coordinate planes are given in 

Figs. 7, 8, and 9. 






Fig. 9. 



Fig. 7. 



Fie. 8. 
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Orbits of Class B. 
Equilibrium point (a). 

e p t = e ->' r y[e+e i ( — 1.20 + 1.20 cos cjt+3. 95 sin ax) + ], 

eq 2 = e-"*y[0.75e— e 2 (0.90 +4. 65 cos en; + 2. 46 sin ctt) + ], 

p 2 =1.833+ ()« 2 + a=1.68. 

TABLE 6. 
€=0.01, y=0.1. 



T 


ea? 2 


eyn 


* 


«7«2 





.0100 





.0010 00 


.000713 


.2 


.0094 


— .0088 


.0007 00 


.0004 77 


.4 


.0078 


— .0166 


.0004 90 


.0003 31 


.6 


.0053 


— .0225 


.0003 42 


.0002 31 


.8 


.0023 


— .0259 


.0002 38 


.000163 


1.0 


— .0011 


— .0264 


.000163 


.000116 


1.2 


— .0043 


— .0240 


.000113 


.0000 820 


1.4 


— .0070 


— .0189 


.0000 776 


.0000 581 


1.6 


— .0090 


— .0117 


.0000 534 


.0000 412 


1.8 


— .0099 


— .0031 


.0000 361 


.0000 289 


2.0 


— .0098 


.0058 


.0000 248 


.0000 203 


2.2 


— .0085 


.0140 


.0000169 


.0000 140 


2.4 


— .0063 


.0207 


.0000117 


.0000 0970 


2.6 


— .0034 


.0250 


.0000 0809 


.0000 0664 


2.8 


— .0001 


.0266 


.0000 0560 


.0000 0452 


3.0 


.0032 


.0252 


.0000 0392 


.0000 0306 


3.2 


.0062 


.0210 


.0000 0275 


.0000 0208 


3.4 


.0084 


.0144 


.0000 0193 


.0000 0141 


3.6 


.0097 


.0062 


.0000 0136 


.0000 0095 


3.8 


.0100 


.0027 


.0000 0094 


.0000 0065 



These orbits are shown in Fig. 10. 




fig. io. 
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Equilibrium point (b). 

ep 2 =e-' , » T y[e+e 8 (0.69— 0.69 cos crx— 2.47 sin err) + ], 

e 22 = e-'* r y[0.40e+ f 2 (0.28+4.60cosffT + 3.07smorT) + ], 

p 2 =3.366+()s 2 +...., «j=2.61. 

TABLE 7. 
f=0.01-, y=0.1. 



T 


tXt 


«J/2 


cyp 2 


«7So 





.0100 





.0010 00 


.0004 49 


.1 


.0097 


— .0103 


.000710 


.0003 70 


.2 


.0087 


— .0199 


.0005 04 


.0002 34 


.3 


.0071 


— .0282 


.0003 59 


.000166 


.4 


.0050 


— .0345 


.0002 55 


.000117 


.5 


.0026 


— .0385 


.000182 


.0000 83 


.6 


.0000 4 


— .0399 


.000131 


.0000 58 


.7 


— .0025 


— .0386 


.0000 933 


.0000 398 


.8 


— .0050 


— .0347 


.0000 670 


.0000 275 


.9 


— .0070 


— .0284 


.0000 481 


.0000190 


1.0 


— .0086 


— .0202 


.0000 345 


.0000131 


1.2 


— .0100 


— .0003 


.0000178 


.0000 063 


1.4 


— .0087 


.0196 


.0000 092 


.0000 031 


1.6 


— .0051 


.0344 


.0000 047 


.0000 016 


1.8 


— .0001 


.0399 


.0000 024 


.0000 0087 


2.0 


+ .0049 


.0348 


.0000 012 


.0000 0047 


2.2 


+ .0086 


.0204 


.0000 006 


.0000 0026 


2.4 


.0100 


.0006 


.0000 003 


.0000 0014 



These orbits are shown in Fig. 11. 




Fig. 11. 
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Equilibrium point (c). 

e p 2 =e-^y[e+e 2 (25— 25 cos err— 23 sin <tt) + ], 

fg, = e-' ,2T y[3.08e+e 2 (77— 56 cos <tt+31 sin <tt) + ], 

p 2 =0.475+()e 2 +...., cr = 1.07. 

TABLE 8. 
e=0.01, y=0.1. 



T 


tx 2 


e#2 


«7P 2 


«7<?2 





.0100 





.0010 


.0032 9 


.1 


.0099 


— .0022 


.0009 3 


.00317 


.2 


.0098 


— .0043 


.0008 7 


.0030 6 


.3 


.0095 


— .0063 


.00081 


.0029 7 


.4 


.0091 


— .0083 


.0007 6 


.0028 7 


.5 


.0086 


— .0103 


.0007 2 


.0027 9 


.6 


.0080 


— .0120 


.0006 8 


.0027 


.7 


.0073 


— .0137 


.0006 5 


.0026 2 


.8 


.0066 


— .0152 


.0006 3 


.0025 4 


.9 


.0057 


— .0165 


.0006 


.0024 6 


1.0 


.0048 


— .0176 


.0005 8 


.0023 9 


1.5 


— .0003 


— .0199 


.0005 4 


.0020 5 


2.0 


— .0054 


— .0169 


.0004 6 


.00171 


2.5 


— .0089 


— .0090 


.0004 2 


.0013 7 


3.0 


— .0100 


.0014 


.0003 7 


.0010 6 


3.5 


— .0082 


.0140 


.0003 


.0007 8 


4.0 


— .0042 


.0183 


.0002 3 


.0005 7 


4.5 


.0010 


.0200 


.00014 


.0004 1 


5.0 


.0060 


.0162 


.00012 


.0003 


5.5 


.0092 


.0078 


.0000 8 


.0002 3 


6.0 


.0099 


— .0027 


.0000 6 


.00019 



These orbits are shown in Fig. 12. 

y 




Fig. 12. 
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